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This manuscript considers multi-criteria based multi-objective approach with 
technical, economic and environmental indices (TEE) for optimal placement 
and sizing of distributed generation (DG) units in the distribution network. 
Technical criteria include indices of active energy losses, voltage deviation; 
whereas economic criteria include the index of cost of DG installation, and 
environmental index considers the various greenhouse gas (GHG) emissions 
from generating unit’s and biomass DG. Combined sensitivity analysis is 
applied for sorting the candidate nodes for DG placement and reducing the 
search space. Multi-criteria decision-making among TEE factors are 
addressed using a scientific approach named Analytic Hierarchy process 
(AHP) approach. The impact of prioritized solutions is analyzed in terms of 
three scenarios formed using AHP in the form of TEE criterion. 
The developed formulation is tested on IEEE 33-bus bus radial distribution 
system and is solved using hybrid optimization approach (hybrid GA-PSO) 
and AHP based scenarios performed better than base case scenario (non- 
prioritized scenario). 
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1. INTRODUCTION 


In the present power scenario, the optimal penetration of distributed generation (DG) units for power 
generation has given multiple benefits in terms of technical economic and environmental (TEE) criterion. The 
most of the centralized power generators utilizes fossil fuels with transportation costs and high the greenhouse 
gas (GHG) emissions. The technical analysis of DG power generation includes power loss reduction, voltage 
profile improvement, and line loading capacity enhancement, network stability improvement and power quality 
improvement. The economic aspect included the cost factors covers the cost of power and energy losses, costs 
related to operation and maintenance of DG and distribution network expansion, DG installations, peak power 
losses etc. However, environmental indices focus on reduction in GHG emissions, penalty due to emissions 


and others keys points. 


Recently, Prasad et al. [1] presented a techno-economic multi-objective whale optimization-based 
approach for placement and sizing of wind and solar based DGs for power loss reduction and voltage profile 
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improvement in the distribution system. Annual economic losses are substantially reduced with reduction in 
annual energy losses. Ali ef al. [2] exercised a genetic algorithm (GA) based approach for DG placement and 
sizing in the distribution network for improvement voltage profile, increasing spinning reserve and reduction 
in line power flows and losses. Raj and Saravanan [3] proposed ‘dwarf mongoose optimization technique’ for 
optimal placement and sizing of DG and STATCOM in the distribution network using loss sensitivity factors 
approach in the radial distribution network. It resulted in reduction of power losses, improvement in voltage 
profile, improvement in system operation cost along with environmental constraint. Sun ef al. [4] developed a 
methodology to examine the impact of DG placement on the efficiency and reliability distribution network of 
with time varying load pattern. Switching operations of DG’s may affect the optimal operation of the 
distribution system. Yaldiz et al. [5] developed a stochastic optimal power flow method for DG integration in 
an active distribution network for evaluation of different indices such as DG integration ratio, generation 
curtailment ratio and adjustable power factor. The impact of DG placement into the distribution system in terms 
of environmental impact assessment has been concluded in [6], [7]. Akella et al. [6] proposed social, economic, 
and environmental effects of renewable energy systems for the grid connected distribution systems. The 
reduction in CO2 emission is determined for both baseline and plants cases is calculated. However, model 
considers only CO2 emissions. Descateaux et al. [7] developed a methodology to assess the implications of 
introducing different carbon tax levels and to assess the greenhouse gas (GHG) abatement performance with 
DG penetration. 

From the above literature, it may is observed that most of the literatures available on optimal sizing 
and placement of DG units considers non-prioritized indices in the formulations. Therefore, it is not possible 
to analyze the impact of different indices in the objective function and relevant results. Therefore, multi- 
objective problems are converted into single objective problem by giving equal weightage to all the objectives. 
A promising logical and scientific technique for weighing factor selection needs further study, as the selection 
of weighing factors decides the quality and optimality of the obtained solution. Further, the reduction of GHG 
emission due to DG placement is also an issue of concern around the globe. This issue has not given due 
consideration in the existing literature. 

Hence, the proposed study aims to determine the optimal placement and sizing of renewable energy 
based DGs considering to technical, economic and environmental benefits by DG. An analytic hierarchy 
process (AHP) technique is used in this multi-objective problem for the selection of weighting factors using 
multiple criteria multi decision process. The weights are mathematically driven and varied by consistency ratio 
test. A hybrid GA-PSO optimization approach is applied to solve the developed formulation. The obtained 
results reveal the effectiveness of the proposed methodology. 


2. PROPOSED METHODOLOGY AND PROBLEM FORMULATION 

The developed mathematical formulation for optimal allocation of dispatchable DG namely micro- 
turbine (MT), gas-turbine (GT), and fuel-cell (FC) based DG types in a radial distribution network is multi- 
objective in nature. Weighting factor are multiplied in accordance to AHP and non-AHP scenarios. The objective 
function comprising of various technical, economic and environmental indices are discussed as follows: 


2.1. Objective function 

The objective function of proposed formulation for optimal allocation of dispatchable DG units in a 
radial distribution network comprises techno-economic criteria along with environmental benefits, which is to 
be minimized. The objective function includes the indices for network losses, voltage deviation, installation 
cost of DG units and GHG emission. The detailed formulation of each term is presented in the subsequent 
sections [8]. 


2.1.1. Index for active energy losses 
Mathematically, the index of active energy losses can be expressed by (1) as [8], [9]: 


ELoss 
IREL = 2258 (1) 


ELoss 


network with DG integration and without DG integration, respectively (in kWh) and computed by (2) and (3). 


[ya Yaar aay Bm (2) 
Bisco OP erp, (3) 
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Where, P47! and P47 are the active power losses corresponding to d‘" demand level before and after DG 
placement, respectively, in the network (in kW); Ng; is the number of load levels; and TDg is the duration of 
occurrence of d‘” load level in a year (in hours). 


2.1.2. Index for voltage deviation 
Mathematically, the index of voltage deviation can be expressed by (4) as [8], [9]. 


w7.B*s{ (v2, -Vmin) + +(VAva- ~Vmax) } 


IVD = 
pole Visi) +(Viva ~Vmax) } 


(4) 


Where, [VD is the index of voltage deviation; Via and V2 ‘q are the voltage magnitudes at if bus corresponding 


to d“" demand level before and after DG placement, respectively, in the network (in kV); and Vinin and Vinax 
are the minimum and maximum permissible limits, respectively, on voltage magnitude of the system (in kV). 


2.1.3. Index for environmental impact 
For better environmental impact, this index should be as minimum as possible. Mathematically, the 
index of environmental impact can be expressed by (5) as [10]. 


EDG 
IE] = -G#¢ (5) 
EGHG 
Where, [EI is the index of environmental impact; and E26, and Egy¢ are the annual GHG emissions with DG 
integration and without DG integration, respectively (in kg) and are computed by (6) and (7). 


ERS. = aL ah] pee. RECO sy ee a DG -{PET - RES + PM? . RE"}]- TD, (6) 
ae = ey Pgh RE® . TDa (7) 


Where, P"? and PS? are the active power imported from the upper grid corresponding to d*” demand level 
before and after DG placement, respectively (in kW); RE®? is the emission rate of the upper prid (in kg/kWh); 


xp @ is a binary variable indicating whether j‘” bus is a candidate location or not; Pat and Pr are the active 


power generation from gas turbine and micro turbine based DG units, respectively, at j*” bus ‘ean eauontine to 
d® demand level (in kW); and RES? and RE™’ are the emission rates of gas turbine and micro turbine based 
DG units, respectively (in kg/kWh). 


2.1.4. Index for installation cost of DG units 
Mathematically, the index of installation cost of DG units can be expressed by (8) as [10]. 


N 
y Bus xPELy ha T+ICMT.p PMT +1CFC.p Fe} 


ie 2 = 


Pa) ok na GS ce ee 8 
ENB PO fic OT BRT, 1M PME, v1CFO PEL] - 

Where, JIC is the index of installation cost of DG units; [C°7, ICM and IC*© are the estan cost of gas 
turbine, micro turbine and fuel cell based DG units, respectively (in sande a 5 pur and PY ‘z are the sizes of 
gas turbine, micro turbine and fuel cell based DG units, respectively, at j* i candidate ioearon (in kW); and 
PST ., PMT. and P£S.. are the maximum sizes of gas turbine, micro turbine and fuel cell based DG units, 


respectively (in kW). Finally, the objective function, C can be expressed by (9). 


min C (9) 
(Pfr PME Pr GPG? of? 2, PhT® al: d=1 to Nat) 
Where: 
C=w, XIREL + w.XIVD + wz XIJEI + wy XIIC (10) 


where, W,, W2, W3, and w, are the weighing factors to the indices for active energy losses, voltage deviation, 
environmental impact and installation cost of DG units, respectively. The weighing factors are selected in such 
a manner so that the (11) is satisfied: 


Wy + W2+ wW3+w,=1 (11) 
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2.2. Constraints 

The objective function with various indices needs to consider constraints too. The developed 
formulation includes constraint such as operating and planning restrictions to be followed during the allocation 
of DG units in a radial distribution network. Both equality and inequality type of bounds considered in the 
optimization are discussed as follows: 


2.2.1. Thermal limit of feeder 

Since each distribution feeder of a distribution network is designed to carry a certain amount of 
maximum current (referred to as ‘Thermal Limit’ of feeder), loading a feeder beyond its thermal limit can 
cause severe damage to it [10] which is governed by (12). 


[(efsy” + (QFS)° < FC; for j = 1to Ngy, — 1, and d = 1to Ny (12) 


Where, Pi * and Qi 5 are the active and reactive power flows, respectively, through sending end of j*" feeder 
corresponding to d‘” demand level after DG placement in the network (in kW and kVAr, respectively); and 
FC; 9 is the apparent power capacity of j th feeder (in kVA). 


2.2.2. System voltage profile 
For a given demand level, the voltage magnitude at every bus should be maintained within its pre- 
specified limits after DG placement as [11] which is governed by (13). 


Vinin S Via S Vinax for i = 2 to Ng,,, and d = 1 to Ng, (13) 


Where, Via is the voltage magnitude at i‘” bus corresponding to d‘” demand level after DG placement in the 
network (in kV). 


2.2.3. Size and operating limit on DGs 
The optimal sizes of different types of DG units are computed by (14), (15), and (16) considering the 
peak load of the system and are restricted to certain fraction of total peak demand of the system as [11], [12]. 


Par < 6° "Qa: yi PP J =1to Neus (id) 
PMT < DC. yous pP j =1t0 Neus ™” 
Pir < xe Qs ee PP J =1to Neus (te) 


Where, a is the fraction of total peak demand of the system. 
For other than peak load operation, the optimal power generation from different types of DG units are 
restricted to their capacities are given by (17), (18), and (19) 


Pep Ss xpos Py j =1to Ng,; andd = 2to Ny (17) 
Pa Sap Py j =1to Ngys and d = 2 to Ny (18) 
Psa e j =1to Ng,; andd = 2to Ny (19) 


2.3. Analytic hierarchy process for selection of weighing factors 

The final objective function given by (10) and (11) is obtained by taking the sum of indices for active 
energy losses, voltage deviation, environmental impact and installation cost of DG units, given by (1) to (10) 
respectively, after multiplying them with suitable weighing factors. Therefore, the selection of proper 
weightage to each index becomes very important to achieve an optimal solution, which results minimum 
investment cost of DG units and maximum technical and environmental benefits. This is because the optimal 
sizes of DG units are computed considering the peak load condition in the network. Mathematically, this is 
expressed by (20). 


W, = W2 = W3 and Wy, W2,W3 > W4 (20) 


Since, the sum of all weighing factors is equal to unity; the following relations can also be given by (21). 
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1-w4 1-w4 
W, = W2 = W3 = and 


> Wy (21) 


The above equation implies to (22). 
W4 < 0.25 (22) 


As the capacity of DG units is obtained for peak load condition, the installation cost of DG units is 
fixed at peak load condition. Therefore, for obtaining optimal DG operation at load levels less than peak load, 
the index for installation cost of DG units is not necessary and is set to zero to (23) as: 


w, = 0 (23) 


the values of weighing factors w,, wz, and wz are computed by using analytic hierarchy process (AHP) 
technique [13], [14]. In the proposed formulation, three technical indices, i.e., real energy losses, voltage 
deviation, environmental impact, in the objective function is taken as the input criterion for decision-making. 
Once the comparison matrix is found to be consistent, the values of weights are computed and used to form 
scenario-2, 3 and 4 reflecting the priority to one out of three technical indices. Scenario-1 is anon-AHP scenario 
and reflects equal importance to all indices [15]. 


2.4. Net annual savings 

The economic analysis of DG integration is calculated in terms of net annual savings (NAS). Where, 
NAS is the net annual saving to the distribution utilities (in $) [10]. It considers annual economic benefit with 
annual investment of distribution utilities. The net annual saving due to integration of DG units in the 
distribution system computed by (24). 


NAS = AB — AI (24) 


2.4.1. Annual economic benefit to distribution utilities 

DG provides active power support to the loads connected at nearby buses, and therefore, reduces the 
energy imported from the upper grid. After integration of DG units into the existing distribution system, the 
annual economic benefit to the distribution utilities due to the reduction of energy imported from the upper 
grid can be computed by (25). 


AB = Cx Yin (P92? — P§?1) - TDg (25) 


Where, AB is the annual economic benefit to the distribution utilities (in $); and Cy is the price of energy 
imported from the upper grid (in $/kWh). 


2.4.2. Annual investment from distribution utilities 

With integration of DG units into the existing distribution system, the annual investment from the 
distribution utilities towards procurement, installation, operation and maintenance of DG units are computed 
by (26). 


ATS aes ICT - Peps r(atryrot Panes ICMT . ate retry? ee 
Jat GtryTeriy 1 (atryTM? 4 


FC 
ec DG -ICFC . prc . retry” 
Ji. (+r TFC 1 


aD tee DG “fy vat {por . Oct + + PMI -OCMT 4 Pe : ocr} : TDa} (26) 


Where, AI is the annual investment from the distribution utilities (in $); r is the interest rate; T°’, T”" and 
T° are the life-times of gas turbine, micro turbine and fuel cell based DG units, respectively, (in year); and 
oc, OC and OC are the operating and maintenance costs of gas turbine, micro turbine and fuel cell 
based DG units, respectively, (in $/kWh). 


2.5. Methodology/solution procedure 

The values of weighing factors for Prioritized scenarios, i.e., W,, W2, and w3 are computed using AHP 
technique. The overall four scenarios are formed including Scenario-1 as non-AHP and other AHP based 
scenarios are numbered 2 to 4 are given in: 
- Scenario-1: Equal weightage to all indices. 
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- Scenario-2: Prioritization to index of active energy losses. 
- Scenario-3: Prioritization to index of voltage deviation. 
- Scenario-4: Prioritization to index of environmental impact. 

It has been observed from the various literature studies that equal weightage to all objectives are 
considered in most of the multi-objective formulations. Hence, scenario-1 has been considered in this work as 
base case scenario. Other scenarios have been also generated giving priority to one index over remaining two 
indices. For scenarios 2 to 4, the weights to the indices are being scientifically calculated using AHP [15]. In 
AHP, the comparison matrix has been computed from the input data from power system experts in terms of 
comparative importance based on respective knowledge and expertise. The comparison matrix is acceptable; 
only if the value of its consistency ratio (CR) is below 0.1 [16], [17]. Table | shows the values of weighting 
factors for scenarios. The solution procedure for getting optimal solution of the proposed formulation using 
Hybrid GA-PSO optimization technique and calculation of weighting factors for scenarios 2 to 4 using AHP 
are shown in Figures |(a) and 1(b). 


Table 1. Values of weighting factors and CR for scenarios 1-4 using AHP 
Scenarios Weighing factors CR. 
Scenario-1 W, = 1/3, w2 = 1/3, w3 = 1/3 - 
Scenario-2.  w, = 0.5934, w, = 0.1284, wz = 0.2764 0.0053 
Scenario-3  w, = 0.1655, w, = 0.6098, wz = 0.2247 0.0904 
Scenario-4 — w, = 0.1488, w, = 0.1603, wz = 0.6908 0.0053 
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¥ 
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scenarios C stop) 


(a) (b) 


Figure 1. Flow chart of proposed methodology (a) steps for AHP calculation and (b) solution procedure using 
hybrid GA-PSO optimization 


3. RESULTS AND DISCUSSION 

The developed methodology is implemented using MATLAB environment and tested on two test 
systems, namely IEEE 33-bus radial distribution systems to determine the optimal sizing and siting of 
dispatchable DG units. For both these test networks, the daily load duration curve is approximated by a 
piecewise linear function. The load levels and corresponding time durations are obtained from [18], [19]. 

In this study, micro-turbine, gas-turbine and fuel-cell based DG units are considered for possible 
integration in both the test distribution networks. The installation and size of different DG are considered from 
[20], [21]. In this study, Price of energy imported from the upper grid is 60 $/MWh, Emission rate of upper 
grid is 632 kg/MWh and interest rate is taken as 10%. For the peak load condition in the network, the weighing 
factors for active energy losses index (wW,), voltage deviation index (wz) and environmental impact index (w3) 
are considered to be equal and higher than the index for installation cost of DG units (w,) satisfying (22). For 
other load levels, w, is set to zero as given in (23) as the optimal sizes of DG units are already computed 
considering the peak load condition in the network. 
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3.1. IEEE 33-bus radial distribution system 

The 12.66 kV, 33-bus radial distribution system is used as test system for implementation of proposed 
topology. The GHG emission for base case is 14.17 kt per annum [7], [22], [23]. To identify the suitable buses 
for DG placement in 33-bus radial distribution network, a sensitivity index is computed using (27). 
Mathematically; the sensitivity of active power losses, P, with respect to active power injection at i” bus can 
be given as [13]. Where, P; denotes the active power injection at i*" bus; AP, is the incremental change in P;; 
and P, (P; + AP;) and P, (P;) are the active power losses with injected power P; and P; + AP;, respectively. At 
each node; DG power injection of 10%, 20% and 30% of total active power demand (TAPD) of the system is 
made and the sensitivity indices are computed using (27). 
PL _ Pi(PitAPi)—PL(Pi) 


ee . a 
Sensitivity of P, = OP; AP; 


(27) 

As seen in the Figure 2 among different buses, buses 18, 17, and 16 have the higher value of sensitivity 
index hence, these buses are identified as the candidate buses for DG placement for placement of micro-turbine 
(MT), gas-turbine (GT), and fuel-cell (FC) based DG units, respectively. The optimal sizes of different types 
of DG units are restricted to 25% of total peak demand of the system (i.e., a = 0.25). This results the maximum 
penetration of each DG type as 900 kW. Considering the peak load condition in the network; firstly, an analysis 
of impact of installation cost of DG units on active energy losses, voltage deviation and environmental impact 
is analyzed. Since the weightage to the index of installation cost of DG units, wy is restricted by (20), its value 
is varied from 0.00 to 0.24 in a step of 0.03 and the impact of w, on other indices considered in the objective 
function are analyzed. For each value of w,, other weighing factors, i.e., W,, Wz, and w3 are computed using 
(21) and optimal sizes of DG units are determined using hybrid GA-PSO based approach. The final solution 
and the values of different indices are presented in Table 1. Since the behavior of different indices is not similar 
with the change in the value of w,, the optimal sizes of DG units corresponding to w, = 0.12 is selected to 
ensure cost-effective penetration of DGs with reasonable technical and environmental benefits [10], [23]. Thus, 
the optimal sizes of micro-turbine, Gas-turbine and fuel-cell based DG units become 0.88 MW, 0.88 MW, and 
0.11 MW, respectively. 

In Table 2, Base case scenario reflects equal priority to all indices in the objective function whereas 
scenarios 2 to 4 are prioritized scenario with respect to the indices formulated in the problem [24], [25]. Further, 
all scenarios result in the positive value of net annual saving, therefore, all scenarios are economically viable 
alternatives for the operation of DG units. Among these four scenarios, scenario-4 results maximum net annual 
saving. Thus, it is beneficial to the system to operate the DG units according to the results of scenario-4. It also 
causes minimum losses in the network, higher economic savings and minimum GHG emission. 


—O— 10% of TAPD —f— 20% of TAPD 30% of [APD 


SENSITIVITY VALUE 
o 9 
— N 
bh 
q 
| 


1 5 9 13 17 21 25 29 33 
BUS NUMBER 


Figure 2. Sensitivity analysis of 33-bus distribution system 


Table 2. Comparison of TEE performance of 33-bus distribution system under different scenarios 
Scenario Active energy losses Active energy from grid GHG emission Net annual saving 


(in MWh) (in GWh) (in kt) (M$) 
Base Case 827.99 22.42 14.17 - 
Scenario-1 522.13 14.89 11.51 17.79 
Scenario-2 523.54 15.10 11.59 17.72 
Scenario-3 562.43 15.63 11.77 15.40 
Scenario-4 487.37 14.54 11.38 19.87 


4. CONCLUSION 

In this study, TEE based analysis for optimal placement and sizing of dispatchable DG units has been 
carried out with mathematical multi-criterion decision making approach for selection of weighting factors i.e. 
AHP. Sensitivity analysis for selection of most sensitive nodes is performed which ultimately reduces the 
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burden of iteration and optimization duration. All AHP based scenarios results in higher saving, better technical 
performance and lowering of the emissions as compared to base case. It gives a better sense to consider 
priority based optimal solutions rather than equal priority indices for DG placements and achieve results that 
are more efficient. 
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